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ABSTRACT 


A study of the vibration of thin laminates in fluid 
is done by the finite element method. The analysis is concerned 
with symmetric orthotropic laminates. FEM is applied to 
both the plate and the fluid domains. The fluid is considered 
to be ideal thus offering no damping on the vibrations . 

The effect of laminate orientations, boundary 
conditions and the aspect ratio on the frequencies of the 
laminates in fluid is studied. Effect of change in fluid 
density on the frequencies is investigated. The in vacuo 
results of the laminates are also presented. 

A study of the 12 d.o.f. and 16 d.o.f. elements 
and the relevant B.Cs. is made and 16 d.o.f. element with 
implied B.Cs. is found to be better. 

A better 13-noded fluid element for gradual coarsening 
of fl ui d mesh is generated by static condensation instead 
of symmetric transformation hitherto used by many 
investigators . This is employed in all the investigations . 



CHAPTER -1 


INTRODUCTION 

Natural frequencies of the underwater structures 
such as submerged pipelines, dams, ship hulls etc. are 
significantly affected by the presence of the fluid around 
them. The inertia of the surrounding fluid reduces the 
natural frequencies of the structuresiT^he density of the 
surrounding fluid is the chief parameter affecting the frequen- 
cies. In addition, the geometry of the structure has an 
important bearing on the vibrational characteristics. Vibration o 
plates has an important application in the ship hull vibration. 
With the advent of composites, these materials are increasing 
used in the marine structures, especially ships and boats . 

Very few studies have been carried out on these structures. 

The present study concentrates on the symmetric orthotropic 
liminates in fluid. The study is limited to thin plates. 

The surrounding fluid is assumed be to incompressible 
and non viscuous. Hence the fluid doesnot contribute to the 
damping of vibrations. 

1 .1 REVIEW OF PREVIOUS WORK : 

A good review of the literature on the structural 
vibration in fluid is given by Muthuveerappan 
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One of the earliest Study on the vibration of plates 
in water was by Lamb /~2_7 which was followed by the work of 
Lewis J~5J who extended the analysis to the ship hull 
vibrations. Lamb studied a clamped circular plate in contact 
with infinite expanse of water on one side by Rayleigh method 
and suggested a formula for fundamental frequency parameter. 

He extended the study to fluid on two sides also. The results were 
proved by Powell through experimental studies. 

Babayev /~5_7 studied the vibration of a cantilever 
plate in water. Lindholm et al £~ 6_7 compared the experimental 
results of the vibration of cantilever plates in fluid with 
the theoretical predictions based on the thin plate theory. 

Clough flJ. Malhotra £~ 8_7, Asghar 9_7 and Fritz /“10_7 
did the analysis through empirical formulae and the conclusions 
of their work were used by J /~1_7 for comparison. 

One of the earliest finite element formulations was 
by Zienkiewicz et al /~11_7 who analysed the vibration of 
complex structures in fluid. Zienkiewicz and Newton /"”12_7 
extended their studies to the coupled vibration with the 
inclusion of the compressibility of the fluid. Holand /~13_7 
extended the analysis to some two dimensional interaction 
problems by developing a few different types of fluid elements . 
Nath ri 4,1 5J performed a finite element analysis of 
vibration of high gravity dams due to the earth quake motion and 
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calculated the hydrodynamic pressure on the dams due to this 
motion. Most of the investigaters mentioned above .idealised the 
fluid as two dimensional region. 

One of the earliest study including the three 
dimensional analysis of the fluid domain was performed by 
Chowdhury jAB/t ij . He employed the finite element method for 
both the plate and fluid domains. He solved for the vibration 
of a steel cantilever plate in water which was assumed to be 
non viscous and incompressible. He compared his results with 
that obtained by strip theory /” 6 _7 and from experimental 
studies. A similar study with compressibility of fluid included 
has been performed by Marcus Z~18_7, Aggarwal et al 1 5_7 
extended the analysis of Chowdhury 6_7 to the simply supported 
plate in water. 

A detailed investigation on the pattern outlined by 
Chowdhury was carried out by Muthuveerappan et al £~ 2@,21_7. 

They have presented the variation of frequency parameter with 
fluid density for various plates /~ 21_7# They later on extended 
the analysis to Skew plates in water /”22,23_7. 

1 .2 OBJECTIVE AND SCOPE OF PRESENT WORK : 

The aim of the present work is to analyse the vibration 
of symmetric composite laminates in fluid by the finite element 
method. The surrounding medium is assumed to be ideal (Nonv'iscous 
and incompressible). Thin and small deformation plate theory is 
employed for the analysis . 
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l\i the second chapter finite element formulation 
for both the plate and fluid region is given. The coupling 
of the plate and fluid is discussed. The generation of the 
added mass matrix is outlined. 

In the third chapter results for the isotropic as 
■well as the symmetric composite laminates vibrating in fluid 
are given* Comparison of the 12 d.o.f. and 16 d.o.f. elements and 
the relavent boundary conditions is done first. The effect 
of the laminate orientations, fluid density and the aspect 
ratio is discus-sed. 

Conclusions are given in last chapter. 


I 
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CH/PTER -2 

FINITE ELEMENT FORMULATION 

This chapter outlines the procedure employed in 
solving the vibration of isotropic and composite plates in 
fluid. Equations of motion of plate are obtained by Hamilton’s 
principle and fluid equilibrium equation by Rayleigh Ritz 
method. The procedure is known in literature £~' 16,21_7, but 
has been included here for the sake of completeness. 

2.1 FEM EQUATIONS FOR PLATE : 

A rectangular plate finite element with the forces and 
moments acting on it is shown in the Fig . 1 . Q x , Qy are the 
shear forces (per unit length) along the y and x -edges of 
the plate respectively, M x , My are the bending moments (per 
unit length) acting on the edges parallel to the y- and x-axes 
respectively ; and M is the twisting moment (per unit length) . 

The Lagrangian L for the element is 

L = T- U+V (2.1 ) 

where 

T = Kinetic energy of the element = JJ 2? p. dx dy 

( 2 . 2 ) 

U = Strain energy of the element = ^ dx dy 

(2.3) 
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V = Work done by applied forces on the element 
= ifqw dx dy + Z” w ( Q x n x +Q y n y) fw x(Mx n x +M xyny) + 


w y (M xy n x +M y n y ) “ 7dB 

in which 

w = Lateral deflection of the plate 
M- = Hass /unit area of the plate 
LMJ = L Mx , LKj= L-w^ - Wyy -2 Wxy j 


(2.4) 


(2.5) 


“s*. 

8M *Y 


3x 

3y 

> 

3M 

3y + 

3M 

_Jg 

3x 


= - ( D i 1 

3^w 

3x^ 

- °12 ^ + 2 
3y^ 

= ~ 2 

3^w 

3x^ 

D 99 + 2D, 

^ 3y^ 


^ - CD, 

*x, - 


3^w 

6 3x 3y 
3^w \ 


a^w . „„ a 2 ^ 

16 ^2 +J 26 ^2 + - d 66 a5Tay 


( 2 . 6 ) 


(2.7) 


where , D.^, ^66 are co_e ^i cien ' ts °f the bending 
stiffness matrix /~D_7 of symmetric orthotropic laminate jT'ZkJJ. 


D 11 °12 D 16 
7_D - 7= D 1 2 D 22 D 26 


D 1 6 D 26 D 66 


( 2 . 8 ) 
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The procedure for obtaining these coefficients is given 
in Section 2.3. For isotropic plate the above matrix 
becomes . 


r°_7= 


Eh* 


12(1 -or) 


1 v 

V 1 

0 0 


0 

0 

1 -v 


(2.9) 


v =Poisson , s ratio 

h = Thickness of the plate 

E = Young’s modulus of the plate material 


Hence , 

L = (Uji rt 2 dx dy - 55 A UU SkI dx dy + 55 < 3. w dx dy + 

A A * 1 J A 


V°y n y ) + w x (M x 1 ^)+ 

B w y (M xy n x +M y V- 7dB (2-10) 

The highest order of derivitive in the double 

integrals in above Eq. .(2.10) is two. Hence compatibility 

upto the first order derivitive is needed. 

The highest order of derivitive in the same Eq.(2.10) 
is three. Hence completeness upto the third order derivitive 
is required. 

The polynomial satisfying the above requirements of 
compatibility and completeness for rectangular element (Fig. 2) 
is bicubic polynomial whose coefficients are evaluated with 
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four degrees of freedom (w, w .w ,w ) at each node, and 
° v 7 x’ y* xy 

is well known /~25, 26_7. 
rC e )_ i M 


w'~' = 1_N J ^w^ 

where |_N J = (_N^ ^ j 

Lwj (ne) = L w 1’ w x1*'y1' w y1’ W xy1’ W 2’ 

Intr oduc ing Eq.(2.11) and 


(2.11) 


(2.12) 


,w. 


xy4 J 


\ M ] -r*j W 

[ k ] 

where /B__7 = - 


(ne) 


(2.13) 


N 1,xx N 2,xx 


•N. 


16, xx 


N i,yy N 2,yy ' 

2N 1,xy 2N 2,xy' 


»N. 


16, yy 


,2N. 


16, xy 


(2.14) 


in Bq.(2.10), Lagrangian becomes. 


l = ^L»j (ne) rn_7 (e) H (n8) - \ Lvj (ne) rxj (e> H Cne) + 

LwJ (ne) ^ (ne W-! (ne) (2.15) 

wherej/~ M7^ = mass matrix of the element = ) N|L-1 J dx dy 

A 

.(e)_ 


Tk_7 


; = stiffness matrix of the element ~ 


[{fvj^LVj /"b_ 7 dx dy 


(2.16) 



9 


^f^( ne )= load matrix of the element = q dx dy 
, , A 

^b^ e - boundary force matrix of the element = 


^ ( (Q x V Q y n y^ ^ N | + K n x + "xy 
B 

+ < M xy VMy V ^> dB 


V \ N ’4 

(2.17) 


Applying Hamiton’s principle, equations of motion are obtained : 

/X7 (e) ^wj (ne) + /“K_7 (e) ^ ne )- ^ ne ) + ^bj (be) 

(2.18) 


The assembled equations for the plate are 


rM_7\w} (n) + TK_7 \v]( n) = (b} (b) (2.19) 

_ __ in) f , ^ ^ are the assembled matrices, 

where fKj, /"M_7, [f] ’\ b j- 

2.2 NOTE ON 12 d.o.f . ELEMENT : 

In addition to the 16 d.o.f, element discussed in 
Section 2.1, the 12 d.o.f. element with w, w and w as nodal 

x y 

degrees of freedom is popularly used in static and dynamic 
analysis of plate. This element does not satisfy the compatibility 
of w along the x-edge and w along the y-edge. For static 

y X 

analysis a comparitive study of these elements is well 
documented, see /"~25,26_7. It is shown that 16 d.o.f. element 
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gives much better results than 12 d.o.f . .element . Comparisons 
between these two elements have been done for equal number 
of elements and not for the comparable size of the solution 
matrix, thus giving the impression that for practically 
acceptable accuracies 12 d.o.f. element is computationally 
cheaper . 

It seems no comparitive study has been done to 
study these two elements for dynamic analysis of plates. 

The present work includes this study. 

2.3 BENDING STIFFNESS MATRIX FOR ORTHOTROPIC LAMINATE : 


For specially orthotropic Lamina, the stiffness 


matrix is 


r^j 


Q. 

Q 

0 

E t 


11 


12 


Q 


12 


Q 


22 


0 


66 


when Q 


11 


Q, 


C 1 -*lt ’W 


V LT E T 


12 


q 22 = 


V TL E L 


( 1 - v t T * tl ) 0 -*lt ^tl } 

S T 


LT TL ; 
E r 


( 1-V LT V Tl) 


Q 66 “ G LT 


( 2 . 20 ) 


( 2.21 ) 


E^, E^, G^rp are the longitudinal, tranverse and shear modules 
of this orthotropic lamina and and are the major and 
minor Poissons ratios. 

Stiffness matrix for the generally orthotropic 
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lamina with fibre orientation at angle 0 to x-axis is 


fSJ- r^J Z~Q_7 7~t_7, 


( 2 . 22 ) 


where 


Cos 2 © 


/“T_7= Sin*© 


Sin 2 © 

Cos 2 © 


Sin© Cos© Sin© Cos© 


2 Sin© Cos© 

2 Sin© Cos© 
Cos 2 © -Sin 2 © 


(2.23) 


Z~Q_7 is changed suitably to connect tensor strains. 

For composite plate of n orthotropic laminates, 

bending stiffness coefficients are, 

D i;f 3 J, h k-b> ^ = 1-2.6 (2.24) 


where h^ is the distance of the upper surface of the k 1 
laminate from the geometric mid surface. 


2.4 FEM EQUATIONS FOR FLUID : 


The differential equation governing the pressure 
distribution during small amplitude motion of an incompressible 
and nonviscous fluid is 


+ ^2 + ix = 0 

3x z 3y^ 3z z 


(2.25) 


where p represents the excess pressure over the hydrostatic 
pressure due to the oscillation of the fluid. 

The corresponding integral to be minimised is /”27 7 


= Z~( H) 2+ ( §y) 2+ (ff) 2 _7 dx d y dz " IS dB 

V B (2 .26) 
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The highest order derivative in the triple integral 
of Eq.(2.26) is one. Hence compatibility of p alone is 
required . 

The highest order derivitive in the Eq.(2.26) is also 
one . Hence completeness upto the first derivitive of p is 
required . 

The polynomial satisfying the above requirements of 
compatibility and completeness for eight noded brick element 
(Fig. 3) is tri-linear polynomial whose coefficients are 
evaluated with only pressure as the nodal degree of freedom 
and is well known in literature r^j 

p (e) = U?J $p} (ne) (2.27) 

LPJ= IP, P 2 P 8 j 

LP) (ne) =|Pi P 2 P 8 i 

Substituting this in the Eq.(2.26), one gets 

I - \ LPj (ne) rH_7 (e) \p^ (ne) - LPJ (ne) ^s} (be) 

(2.28) 

where /“h_7 ^ = Fluid stiffness matrix of the element = 

\ p >xlLP. x J ^P.y^ ^P.y J+ [P. z ^LP, x J)dx dy dz 
V 


(2.29) 
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Boundary matrix of the fluid element = 


H dB 

B 


(2.30) 


Using Rayleigh Ritz method, Eq.(2.28) gives the 
fluid equilibrium equations 

ZX 7 Ce) \p} (ne) = [g]^ (2.3 

Assembling the element matricies 


tnj 

where £~H_7 , (gl ^ '' are assembled matrices. 


(2.32) 


In addition to the 8-noded fluid element mentioned 
earlier, 13-noded fluid element with pressure as the d.o.f. 
is employed for gradual coarsening of the fluid finite element 
mesh. The procedure for deriving the stiffness matrix for this 
element is discussed in Section (2.5). 


2.5 13-NODED FLUID STIFFNESS MATRIX : 


A 13-noded fluid element is shown in Fig. 4. This 
element has nine nodes on one surface and four on every other 
surface. To derive the stiffness matrix for this element, 
eight of the 8-noded elements are first assembled (Fig. 5) 
to form a (27 x 27 ) matrix /~H_7 Chowdhury /”16_7 employed 
a transformation matrix which preserves the linear variation of 
pressure in the z -direction for eleminating the 14 of the 
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unwanted degrees of freedom. For more accurate results 
static condensation technique has been adopted in this work 
to get the required stiffness matrix from the assembled (27 x 27 ) 
/~H_7 matrix. The procedure is outlined below. 


The assembled equilibrium equations for eight, 
8-noded fluid elements are 


H 


pv (ne) - <s 


(be) 




Partitioning the above matrix 

(ne) , „ N (be) 


H 11 H 1 2 


H 21 H 22 



Si 

. S 2 


where ^p.^ = required 13 d.o.f. 

\]? 2 ^ = Unwanted 14 d.o.f. 


(2.33) 


(2.34) 


Hence /"~H 

H (ne) +&H 2-7 ^ (na) - 


r h 21 _7 

^4 (ne) = 

<M 


(be) 


(2.35) 


(be) 


(2.36) 


From Eq. (2.36) 


{v z ] ne) = - r« 22 _7 - 1 rH 21 _7 ^V“ >+ fH 22 _7 - 1 [g 2 ^° 


(2.37) 
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Substituting in Eq. (2.35) 

( p*^j- r^ 2 j rx zz j '' 1 / tj 21 _7 (ne) - 

(g^ (be) - Th 12 _7 Z~ h 22 _7 - 1 ^ (be) 

(2.38) 

Thus stiffness matrix for 13-noded element is 

rn_7 («) = ra^j- p» zz j - 1 rn 21 _/ 

(2.39) 


2.6 COUPLING OF PLATE AMD FLUID EQUATIONS : 


The equation of motion for the plate Eq. (2.19) is 

raj r^j W (n) • ^ (n> + M (n) 

(2.40) 


From Eq. (2.17) 

^f^ ne ) = \ q dx dy 

A 

The pressure p of the fluid acts as a distributed force q, thus 

£f^( ne )= $^N$ ^ dx dy (2.41) 

A 

Substituting jf e ^ from Eq. (2.27), one gets 

^f^ ne )= LpJ dx dy ^p^ ne )= /"L_7^ ^p ^ 

A 


(2.42) 
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where 


r Lj7 (e) = Interaction matrix for the element = 

^ LpJ dx (2.43) 

A 

Assembling the ^f^( ne ' S> one gets 

l (n) = Z“l_7 \p^ (n) (2.44) 

where /~L_7 is the overall interaction matrix, and ^pj^ n ^ 
is the matrix of pressures at the interacting nodes. 

Substituting for ^p^ n ^ from Eq. (2.32) in Eq.(2.44), one gets 

^ (n) = Tl_7 HiJ - 1 <»g] (b) (2.45) 


where J~ H_7 _1 is the fluid stiffness matrix inverted after 
applying the boundary conditions # It may be noted that after 
inversion of £^~ H_J7 matrix, only part corresponding to the 
interaction nodes is used in Eq. (2*45) . 

Matrix is given by Eq. (2.30) and is 

^gf e) - \ IgciB (2.46) 

B 

At the solid fluid interface 

|E = - p, w (2.W) 

8n r f 

where p^? is the density of the fluid and w is the plate 


acceleration. 
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Using Eq.(2.46) and Eq.(2.47), one gets 

H (be) " -p f «' (e)dB 

B 

Using Eq.(2.11), this becomes 

^ < be )= - p f ^ P | L N J dB (ne) 

= - p f /“k_7 {e)T (ne) 

Assembling over the plate, one gets 

\s] w = - p , (n) 

Substituting this in Eq.(2.45), one gets 

^ (n) = p f Tl_ 7 ZX7 _1 r L _7 1 [«} (n) 

Thus Eq.(2.4o) becomes 

( ni_7+ yT M A-7)H (n) + ZX7 H (n) = 1 


where is the added mass matrix 

D\J - Pf Z~ L _7 rn_7 ■ 1 ZX7 T 


(2.48) 


(2.49) 


(2.50) 


(2.51) 


b^ n) 

(2.52) 


(2.53) 
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CHAPTER -3 

RESULTS AND DISCUSSIONS 


Firstly ,a comparison of 16. d.o.f. and 12 d.o.f. 
elements along with relevent boundary conditions is made. 

This is followed by a study of isotropic plate in fluid. A 
detailed investigation of the symmetic composite laminates in 
fluid is presented thereafter. 

3.1 PLATE ELEMENTS AND BOUNDARY CONDITIONS : 

Frequencies for the simply supported plate employing 
16 d.o.f. and 12 d.o.f. elements are given in Table 1 and 2 
for classical as well as classical with implied B.Cs. Similar 
results for the cantilever plate are given in Tables 3 and 4, 
and for clamped -free -clamped free plate in Tables 5 and 6. 

Table 7 shows comparison of both elements for solution matrices 
of comparable sizes. 

Before discussing these elements further, the B.Cs. are 
catalogued. The classical B.Cs. for various edges one : 

Clamped edge : For edge x=a to be clamped 

w I = 0 ; f£ I =0 (3.1) 

1 x=a ' x=a 


Simply supported edge : For edge x=a to be simply supported 


=0 



=0 


w 


x=a 


f 


x=a 


(3.2) 
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Free edge : For edge x=a to be free 




x 


x=a 


3M xv 

”°* (Q x + W 21 - ) 


=0 


x=a 


(3.3) 


The implied B.Cs are 
Clamped edge : For x=a to be clamped 


3w 

3y 


= 0 , 

x=a 


3^w 
3x §y 


x-a 


0 


(3.4) 


Simply Supported edge 
supported 


3w 

3y 


» 0 

x=a 


For edge x=a to be simply 


(3.5) 


The frequencies for isotropic plate are presented 
in the form of non-dimensionalised frequency parameter X, 

\=ua 2 J^ (3.6) 

Where w = circular frequency 
a = length of the plate 

it = mass per unit area of the plate = ph 
p = density of plate 

D = flexural rigidity of plate = — — 

12(1 -v^) 

h = thickness of plate 
v = Poisson’s ratio. 

Tables 2,4 and 6 show the results for the 12 d.o.f. 
element for different B.Cs. It is observed that the frequencies 
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by classical B.Cs. and classical with implied B.Cs. are 
essentially same. But in both these B.Cs., convergence is 
monotonically increasing in some modes, is monotonically 
decreasing in other modes; and sometimes the convergence is 
non-monotonic . For the simply supported case for which exact 
results are known it is seen that the FEM frequencies are 
less than the exact results . 

Tables 1,3,5 show the results for the 16 d.o.f. 
element for classical and classical with implied B.Cs. 
Comparison of the results show that the frequencies are very 
accurate by classical with implied B.Cs. Sometimes the 
results with classical B.Cs. only show an appreciable error. 
The results for the former B.Cs. are generally found to be 
more accurate even with a fewer elements. The convergence is 
always monotonically decreasing, where as with classical B.Cs. 
only, convergence is sometimes monotonically increasing and 
some times monotonically decreasing. It may be noted that the 
trend of monotonic decrease in convergence is consistent with 
consistant FEM formulation. 

From the observations above it is clear that 
classical with implied B.Cs gives greater accuracy for the 
16 d.o.f. element. In addition introduction of implied B.Cs. 
reduce the size of the solution matrices thus decreasing the 
computational effort. Hence implied B.Cs. are always used 
in the subsequent discussion and results in this work. 
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Tables 1 and 2 show the results for the two elements 
for the simply supported plate for which exact solution is 
available. It is evident that the 16 d.o.f. element results 
are better than that of 12 d.o.f. elements. Even with a fewer 
elements, the 16 d.o.f. element gives quite accurate results. 
The discrepency between the two elements is appreciable in 
the fourth mode where the 12 d.o.f. result shows an error of 
9 percent where as for 16 d.o.f. it is almost zero. With the 
16 d.o.f. element the convergence is always monotonically 
decreasing where as for the 12 d.o.f., element convergence is 
monotonically increasing in some modes, monotonically 
decreasing in some other modes and in some modes not monotonic 
at all. 

Tables 3 and 4 give similar results for cantilever 
plates . Same trends as those for simply supported plates are 
observed. But the results by 12 d.o.f. element are also quite 
accurate, thus showing that the clamped edges are less affected 
by the slope incompatibility of the 12 d.o.f. element than 
the simply supported edges . 

Tables 5 and 6 show similar studies for the clamped - 
free -clamped -free plate, and same trends as those for the 
cantilever plate are observed. 

Table 7 shows the comparison of the two elements for 
the comparable size of solution matrices * It is obvious that 
the 16 d.o.f. element is superior to 12 d.o.f. element 
computationally . 
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From the foregoing discussion it can be concluded 
that the 16 d.o.f. element with implied B.Cs. is more accurate, 
computationally more efficient and gives monotonically 
decreasing convergence. Hence this element is employed in 
the further investigation in this work. 

3.2 COMPUTATIONAL DETAILS : 

The spread of fluid around the plate and the depth 
below it are in reality infinite in extension. But the effect 
of fluid is negligible after a certain distance from the plate. 
Hence it is sufficient to take a finite extension of the fluid 
for analysis. Muthuveerappan /~1_7 has shown that it is 
sufficient to take the fluid spread around the plate to be 
equal to the maximum plate dimension and the water depth below 
to be 2.75 times the maximum plate dimension. 

The fluid domain with the finite element discritisation 
is shown in the Fig. 6. The layer by layer description of these 
fluid elements is shown in Fig. 7. In general 8-noded tri -linear 
elements have been used. For gradual coarsening of the mesh the 
13-noded element derived by static condensation has been 
employed. Eighty eight fluid elements with 186 nodes including 
51 surface nodes have been used. Out of these 88 elements, 

16 elements are 13-noded. 

The overall interaction matrix coupling the plate with 
the fluid degrees of freedom is obtain by assembling the 
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elementary interaction matrices of Eq. (2.43). In the 
presence of fluid on both sides, this matrix is evaluated 
considering the interfaces above and below the plate. The 
added mass matrix accounting for the inertia of the fluid 
is evaluated from the interaction matrix and the fluid 
stiffness matrix from Eq.(2.53). This is then added to the 
plate mass matrix and the resulting Eq.(2.52) is solved for the 
eigen values after applying the B .Cs . A NAG F02BJF routine is 
employed for eigen value calculations . 

3.3 ISOTROPIC PLATE IN FLUID : 

Tables 8 and 9 give the results of the present work 
for isotropic cantilever and simply supported steel plates 
vibrating in water along with the earlier results. The plate is 
at a depth of 0.25 times the length of plate. The results of 
the present work include both types of plate elements 
previously discussed in Section 3*1 . Results employing the 
1 3-noded fluid element derived by static condensation are 
compared with the results by 1 3-noded element used by /“16,2lJ/. 
In vacuo results are also presented for comparison. 

Table 8 gives the results for the cantilever plate . 

To check the present computer programs, first the results are 
computed for 12 d.o.f. element with 13— noded fluid element used 
by /"1 6, 21__7. A good agreement is found with these references. 
Also results are computed for this element with the more 
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accurate /“16_7 present 13-noded fluid element. The frequencies 
are found to be slightly less than the earlier ones and 
results are found to be closer to the experimental value. Next, 
results are computed for 16 d.o.f. element with the latter 
fluid element, and results indicate the correctness of the 
present computer programes for this element. 

The simply supported plate results are presented in 
Table 9. The results from the present work for both the types of 
13-noded elements are compared with previous results /~1_7 
obtained by 12 d.o.f. element and 13 noded fluid element of 
£“' 1 6 , 21 J , It is observed that higher mode frequencies are 
slightly higher than those of the 1 _7 . This difference is quite 
appreciable for the fourth mode. These trends are consistent 

f 

with the in vacuo case of simply supported plates . Hence the 
present frequencies are more accurate. 

From the discussion above it can be said that 16 d.o.f. 
element with 13-noded fluid element obtained by condensation 
gives better results. This is adopted in future investigations. 

3.4 COMPOSITE SYMMETRIC LAMINATE RESULTS : 

For composite laminates the frequencies are presented 
in the form of non -dimens ionalised frequency parameter, X, 

X = cca 2, (3.7) 

] E T n 

where E^ is the transverse modules of the lamina. 

3.4.1 Composite Plate In Vacuo Frequencies : 

Table 10 lists the in vacuo composite results for 
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a single 0°-lamina and a O/90/O laminate. They have been 
compared with exact results, Jones /“24_7. 

The results for both lamina and the laminate are very 
close to the exact value even upto the fifth mode. The results 
obtained are slightly higher than the exact values as expected. 

It is clear from the above observations that the 
present programme is accurate for the analysis of symmetric 
composite laminates. 

3.4.2 Composite Laminates in Fluid : 

Tables 11-13 give the in vacuo and in fluid frequencies 
of symmetric composite laminates for different boundary 
conditions and various laminate orientations . The laminates are 
at a depth of 0.25 times the length. A typical graphite -epoxy 
composite has been considered for the analysis . The properties 
of this are given in Table 10. 

It is observed from tables that the first mode 
frequency decreases by 55 */. and second mode by 45 */ for all the 
laminate orientations for the three boundary conditions studied. 
For higher modes the decrease depends on the laminate orientations 
and boundary conditions. 

3.4.3 Variation in Frequencies with change in Fluid Density : 

Table 14 shows the effect of fluid density on frequencies 
for the square contilever laminate under investigation. 
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As expected, the effect of the fluid decreases with 
decrease in fluid density . This is due to the fact that the 
inertia of the fluid on the plate decreases with decreasing 
density of the fluid thus lessening the added mass on the 
plate. From the table it is clear that the lower modes are 
affected to a greater extent than the higher modes. The first 
mode is the most affected with a reduction in frequency of 60 '/. 
for plate to fluid density ratio of one~ which gradually 
decreases for higher ratios and is 24 */. for a ratio of 8. 
Similar trend is observed for higher modes . 

3.4.4 Variation in Frequencies with Aspect Ratios (b/a) : 

Table 15 shows the effect of the aspect ratio on the 
frequencies of a simply supported graphite epoxy laminate . 

From the table it is obvious that with decrease in the aspect 
ratio the effect of fluid on the frequencies increases for 
the first two and fourth modes . No particular trend is noticed 
for the other modes. 



27 


VO 

r- 

t*- 

in 

m 

T~ 

r- 

m 

CO 

r- 


4 - 

-4* 

cn 

00 

00 


o 

T~ 

rn 

T- 

*4- 

4 

n- 


4 

00 

CO 

VO 

m 

CM 

VO 

VO 

CO 

m 

VO 

in 

UT\ 

r- 

o 

• 

• 

• 

* 

• 

Ch 

00 

CO 

T~ 

r- 


*4 

•4" 

o 



v~ 

co 

CO 

VO 

r- 


CO 

00 

in 

co- 

r- 

V 

V“ 

m 

o 

o\ 

CTv 

cn 

i> 

oo 

X 

4 - 

4 “ 

r- 

Cn 


4 

in 

in 

VO 

o 

n- 

m 

m 

cn 

r- 

« 

« 

• 

• 

• 

CT\ 


G\ 

oo 

oo 

v— 

4- 

4- 

n- 

<n 


CM 

VO 

VO 

i> 

vo 

4 

n- 

r- 

CO 

00 

r- 

4 

4* 

r— 

v~ 

• 

• 

♦ 

• 

• 

cn 

cn 


c n 

O 

v~ 

4* 

4- 

i> 

o 

v- 


00 

r- 

c- 

VO 

VO 

4" 

m 

m 

00 

o 

n- 

n- 

n* 

4 

VO 

• 

* 

• 

• 

* 

CTn 

cn 

CT\ 

cn 

r- 

v~ 

4- 

4* 

n* 

o 

V* 


C\J 

v~ 

v- 

in 

V 

00 

cn 

cn 

O 

O 

n- 

o\ 

cn 

<T\ 

n- 

• 

« 

* 

* 

■ 

<n 

CM 

cM 

m 

co 

v* 

in 

in 

CO 

V“ 

r— 


v - 

CM 

m 

4” 

in 

<4 

,4 

4 

4 

4 



Table 2 : Simply Supported Plate - In Vacuo Frequencies - 12 d.o.f. Element 
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Table 5: Clamped -Free -Clamped -Free Plate-In Vacuo Frequencies- 16 d.o.f. Element 
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Table 6 : Clamped -Free -Clamped -Free Plate - In Vacuo Frequencies -12 d.o.f. Element 
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Table 7 : Comparison of In Vacuo Frequencies by two 

Elements for Comparable the Solution Matrices. 
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Table 8: Comparison of Frequencies of a Square Cantilever Isotropic Plate 
in water. p=7861 .11 kg/ m 3, p^=1000 kg/iP, (h/b)=0 .01 31 .All 
results for 16 Elements. 
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Table 9 : Comparison of Frequencies of Square Simply 
Supported Isotropic Plate. Data as in 
Table 8. 



Muthuve erappan 

Present work 

16 d.o.f. 

In Vacuo 
Exact /28 _/ 

X 

, 22 j 

12 d.o.f . , 

/T6_7l 3-noded 
fluid Element 

Present 13- 
noded Fluid 
Element 

N 

7.86 

8.039 

7.856 

19.74 

X 2 

23.92 

24.949 

24.495 

49.35 

X 3 

24.37 

25.518 

24.916 

49.35 

H 

41 .23 . 

45.554 

44.786 

78.96 

x 5 

- 

60.645 

59.752 

98.70 
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Table 10 : In Vacuo Frequencies of a Square Simply 
Supported Composite Plate. %/=20GPa, 
Eij,=2.GPa , Gjjj, = 0 • 7GPa t ' 



Single 0° ! 

Lamina 

0/90/0 Laminate 

X - 

FEM 

Exact /?A_7 

FEM 

Exact /~24_7 

N 

10.371 

10.369 

10.371 

10.369 

X 2 

16 .830 

16.797 

18.019 

17.979 

X 3 

30 . 446 

30.032 

33.975 

33.482 

x 4 

37.443 

37.304 

36.886 

36.749 

Xr 

41 .616 

41 .476 

41 .616 

41 .476 


Table 11: Frequencies, of Square Cantilever Graphite Epoxy Laminate. 
E l =159 GPa, E, r =10.9 GPa, G LT =6.4 GPa, v LT = 0.38, 
a =4m, h =0.3m, p=16l0 kg/ id ?, p^-1000 kg/m? 



'Table 12 ; frequencies of Square Simply Supported Graphite Epoxy 
Laminate Data as in Table 11 



Table. 13 : Frequencies of Square Simply Supported -Clamped -Free -Clamped 
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Table 14 : Effect of F] uid density on Frequencies of 
Square Cantilever Graphite epoxy Laminate 
(45/0/45) in Fluid. Data as in Table 11 . 


<H 

CL 

1 .0 

2.0 

4.0 

5.0 

8.0 

Oo 

(In Vacuc 

N 

0.682 

0.901 

1 .138 

1 .213 

1 .358 

1 .788 

X 2 

2.224 

2.858 

3.481 

3.663 

4.001 

4.867 


5.285 

6 .706 

8.042 

8.419 

9.099 

10.720 


6.628 

8.221 

9.639 

10.028 

10.722 

12.397 

^5 

10.643 

13.167 

15.387 

15.990 

17.656 

19.538 
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Table 15 


Frequencies of Rectangular Simply Supported 
Graphite Epoxy Laminate (45/0/45j .E, =1 59GPa, 


E,p=l0,9 GPa, 
p=2000 kg/m 5 , 


G lt = 6.4 GPa, v lt = 0.38, 
p £ =1000 kg/m’ and (h/b)=0.075. 



In 

Vacuo 


In Fluid 

1 .0 

0.75 

0.3 

1 .0 

0.75 

o.: 

N 

12.931 

17.668 

30.913 

6.070 

7.698 

12.: 

X 2 

24.927 

32.277 

47.403 

14.503 

1 7.634 

23.-1 

X 3 

37.039 

52.322 

72.011 

21 .244 

29.742 

42. 4 

x 4 

40 .01 4 

52.726 

102.308 

26.382 

33.309 

59.6 

X 5 

59 .455 

76.624 

111 .640 

41 .695 

49.630 

76.8 
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No. of plate bending elements = 16 
No. of 8-noded fluid elements = 72 
No. of 13-noded fluid elements = 16 
No. of structural d.o.f. =100 

No. of fluid d.o.f. =186 


Fig. 6 Finite element idealisation of plate and fluid . 
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6 13-r>oded elements 
46 8- noded elements 



4 13- noded elements 
20 8 -noded elements 



6 13 -noded elements 


6 8 -noded elements 



Fig. 7 Layer by layer description 


of the fluid domain. 
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CHAPTER- 4 
CONCLUSIONS 

Based on the results obtained in the previous 

Chapter the following conclusions are drawn. 

(1) The 16 d.o.f. element with implied B.Cs. also is 
more accurate, computationally more efficient and 
gives monotonically decreasing convergence . 

(2) Condensed 1 3 -nod ed .fluid element gives greater accuracy 
than the earlier element ri6_7. 

(3) The effect of fluid is maximum in the first mode 
in all laminates studied . 

(4) For the first two modes, laminate orientations and 
the three B.Cs. studied have practically no effect 
on the percentage decrease in frequencies due to the 
fluid. The effect on the higher modes due to the fluid 
show a dependence on the laminate orientations and 
B.Cs . 

(5) Effect of the fluid, as expected, decreases with 
decrease in fluid density. 

(6) The decrease in aspect ratio (b/a) increases the effect 
of fluid on the frequencies in the first two and fourth 
modes for the laminate studied. Similar trend is not 
observed in other two modes. 
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